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Abstract 

Diffractive photoproduction of vector mesons, 7£> — > VY, where Y is a proton- 
dissociative system, has been measured in e + p interactions with the ZEUS de- 
tector at HERA using an integrated luminosity of 25 pb -1 . The differential cross 
section, da/dt, is presented for —t < 12GeV 2 , where t is the square of the 
four-momentum transferred to the vector meson. The data span the range in 
photon-proton centre-of-mass energy, W, from 80 GeV to 120 GeV. The t distri- 
butions are well fit by a power law, da/dt oc (— t)~ n . The slope of the Pomeron 
trajectory, measured from the W dependence of the p° and <fi cross sections in 
bins of t, is consistent with zero. The ratios dcr 7P _>0y/dt to dcr yp ^ p o Y /dt and 
d(T 7P _^j/^y/dt to da^p^poy/dt increase with increasing —t. Decay-angle analyses 
for p°, <p and J/i/j mesons have been carried out. For the p° and <fi mesons, 
contributions from single and double helicity flip are observed. The results are 
compared to expectations of theoretical models. 
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1 Introduction 



Studies of the elastic production of vector mesons in electron-proton interactions [1] show 
that at large Q 2 , the exchanged-photon virtuality, or at high vector- meson mass, the cross 
section increases with energy faster than is observed in hadron-hadron interactions [2,3]. 
The latter increase is well described by Regge theory [4]. Although the energy dependence 
of vector-meson production in ep interactions can be described by models based on Regge 
phenomenology [5] and the vector dominance model (VDM) [6], it can also be explained 
by models based on perturbative QCD (pQCD) [7,8]. In this case, it is related to the rise 
of the gluon density in the proton as x decreases, where x is the Bjorken scaling variable. 
In pQCD models, a perturbative (hard) scale can be provided by either a high Q 2 or a 
large meson mass. It is also predicted [9] that, in diffractive dissociation of hadrons, the 
squared four-momentum transfer, t, may serve as a hard scale. 

The subject of this paper is vector-meson (V) photoproduction at high —t, which is dom- 
inated by the proton-dissociative reaction, jp — > VY, where Y is the dissociated hadronic 
system. The data cover the photon-proton centre-of-mass energy range 80 < W < 120 
GeV and extend from -t = 1.1 GeV 2 up to -t = 12 GeV 2 for the p°, to 7GeV 2 for the (p 
and to 6.5 GeV 2 for the J/ip meson. The vector mesons were identified via their decays to 
two oppositely charged particles: p° — > tt + tt~, <p — > K + K~ , and J ' /ip — > Z + /~(e + e~, 



2 Theoretical approaches 
2.1 The Regge model 

In Regge theory, diffractive processes are assumed to proceed through the exchange of the 
Pomeron trajectory. The vector dominance model, schematically indicated in Fig. 0(a), 
together with Regge theory, gives the following form for the double-differential cross 
section for the reaction ^p — > VY: 

dM 2 dt ~ nt) W 2 \M 2 ) {My) ' {1} 

where My is the mass of the diffractively produced hadronic state, Y. 

Assuming a linear form for the Pomeron trajectory, a(t) = a(0) + a't, fits [10, 11] to 
the elastic photoproduction of p°, <p an d J /if) mesons in the range —t < 1.5 GeV 2 gave 
a(0) = 1.096 ± 0.021 and a' = 0.125 ± 0.038 GeV" 2 for the p°, a(0) = 1.081 ± 0.010 and 
a' = 0.158 ± 0.028 GeV" 2 for the 0, and a(0) = 1.200 ± 0.009^;^ and a' = 0.115 ± 
0.018to!oi5 GeV" 2 for the J/ip meson. These values may be compared to a(0) = 1.08 [2], 
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a(0) = 1.096 [3] and a' = 0.25 GeV -2 [2], obtained from fits to hadron-scattering data at 
—t < 0.5 GeV 2 . These numbers indicate that the Pomeron trajectory is not universal. In 
fact, it has been suggested [12] that, in the production of vector mesons, an anomalous 
Regge trajectory exchange, with a(0) ~ 1 and a' ~ 0, gives a dominant contribution to 
the differential cross section at large —t. 



2.2 pQCD models 

In models based on pQCD [13-18], the reaction 7p — » VY is viewed as a sequence of 
three successive processes, illustrated in Figs. |l](b) and 0(c): the photon fluctuates into a 
qq pair; the qq pair scatters off a single parton in the proton by the exchange of a colour 
singlet; the scattered qq pair becomes a vector meson and the struck parton and the proton 
remnant fragment into a system of hadrons. The probability of the photon fluctuating 
into a qq pair is parameterised by the photon wave-function. Many models assume that 
the interaction of the qq pair with a parton in the proton is mediated in the lowest order 
by the exchange of two gluons [13, 15, 17, 18]. The exchange of the gluon ladder has 
also been computed [13,14,16,18] in the leading logarithm approximation (LLA). The 
transition of a qq pair into a meson is, however, a non-perturbative phenomenon that 
must be parameterised by the meson wave-function derived from lattice calculations and 
sum rules. 



In the perturbative regime, the cross section for the process 7p — > VY is given by 

= g{X) _ t) ^pL + J2(q(x, -t) + q(x, (2) 



where x is the fraction of the proton momentum carried by the struck parton and g(x, —t), 
q(x, —t) and q(x, —t) are the gluon, quark and antiquark density functions, respectively. 
The cross section for the production of a vector meson through the interaction with a 
gluon (Fig. 0(c)) is about five times larger than through the interaction with a quark 
(Fig. 0(b)). At low x, therefore, the process 7p — > VY should be dominated by the 
reaction in which a gluon from the proton is struck. 

In the two-gluon approximation, the cross section of Eq. (fj) does not depend on W for 
any given t, i.e. a' = 0. In the LLA, this cross section has a power-law dependence on 
W, although a' is expected to be small, a' < 0.1 GeV" 2 for — t > 1 GeV 2 [16], decreasing 
with increasing —t and reaching zero at asymptotically large —t values. 

A previous ZEUS publication [10] showed that da/dt for the 7p — > p°Y and 7j> — > <pY 
processes were larger than those obtained by perturbative calculations performed in the 
two-gluon approximation [15]. These calculations have been further developed by Ivanov 
et al. [17], who proposed that a large contribution would arise from qq fluctuations in a 
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chiral-odd spin configuration. In this case, the quark mass appearing in the coupling 
is interpreted as a current quark mass, and therefore the quark mass in the light-meson 
wave-function is neglected. In such an approximation, da/dt for p° and </> mesons should 
have the same shape for — t > 1 GeV 2 . In addition, the cross sections should exhibit a 
power-law dependence on t, da/dt oc (— t)~ n ; for x > 0.01, the exponent n « 3.8 can be 
estimated from the Ivanov et al. calculation [17] for a fixed value of the strong coupling 
constant, a s , or n ~ 4.8 using the leading-order form for a s (—t) [19]. 

In another two-gluon calculation, that of Forshaw and Ryskin [13], the quark mass ap- 
pearing in the coupling is interpreted as the constituent quark mass. Such models 
predict a dip in da/dt at —t ~ My-, where My is the meson mass. 

A LLA BFKL [20] calculation has been carried out by Bartels et al. [16] for heavy mesons. 
The results are expanded in terms of In (xW 2 /Wq), where W$ = M v — t is assumed. This 
model predicts that da/dt oc (— t)~ n , where n is a function of the ratio —t/M v . For 
the t range covered by the J ftp analysis of this paper, an average value of n ~ 1.8 is 
predicted for a fixed value of a s . The previous ZEUS measurement of the 7p — > J/ipY 
cross section [10] shows that the prediction of this model, using a s = 0.2, agrees well with 
the data. However, the uncertainties due to the choice of a s and the Wo scale are large. 

In a recent paper, Forshaw and Poludniowski [18] have extended the BFKL models of 
Forshaw and Ryskin [13] and Bartels et al. [16], and have given cross sections for the p°, 
4> and J/ip mesons. In this model, which uses a fixed value for a s , a delta-function meson 
wave-function was used and the quark mass appearing in the 'jV coupling was assumed 
to be the constituent quark mass for all three mesons. 

2.3 Helicity structure and cross-section ratios 

The hypothesis of s-channel helicity conservation (SCHC) [21] for vector-meson produc- 
tion states that the helicity of the final-state meson will be equal to that of the initial 
photon. Studies of the elastic photoproduction of light vector mesons at low —t [1] show 
that SCHC holds in these soft diflractive processes. 

In VDM, the cross section for vector-meson production is proportional to the product of 
the coupling of the photon to the vector meson and the meson-proton scattering cross 
section. The 7 — > V coupling depends on the meson leptonic width, T v ^ e + e -, which 
depends on the meson mass, the quark assignment and the wave-function. The SU(4) 
prediction, which ignores the differences in the masses and the wave-functions of the 
mesons, is that the ratio of the production cross sections for vector mesons will be 

p° : uj : <p : J/xjj = 1 : 1/9 : 2/9 : 8/9. (3) 
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In pQCD, the helicity of the final-state meson and the cross-section ratios are sensitive 
to the photon polarisation and the meson wave-function. In the trivial case of a meson 
wave-function given by a delta function, which is a good approximation for heavy vector 
mesons, the helicity of the final-state vector meson is equal to that of the initial photon, 
and thus SCHC will hold. In this case, the cross-section for transverse photons is given by 
da/dt oc MyT v ^ e + e ~ [13,16,18]. Based on the measured values of My and T v ^ e + e - [22], 
the ratio of the photoproduction cross sections for various vector mesons at asymptotic 
values of —t becomes 

p° : uj : <p : J/ if) = 1 : (0.8 x 1/9) : (2.1 x 2/9) : (56 x 8/9), (4) 

which is significantly higher than the SU(4) prediction of Eq. ([3]) for the <fi and J / ip meson. 

A wave-function more appropriate for the light mesons [15, 17] leads to the production 
of mesons in the helicity-0 state by transverse photons and therefore to a violation of 
SCHC. The ratios of the production cross sections for mesons in the helicity-±l states 
by transverse photons are equal to those of Eq. The production of mesons in the 
helicity-0 state by transverse photons, which is expected to be dominant at large —t, is 
given by da/dt oc MyTy^ e + e - rather than M v Ty^ e + e - mentioned earlier. In this case, 
the ratios of the photoproduction cross sections at large —t should be 

p° : uj : : J/tp = 1 : (0.8 x 1/9) : (1.2 x 2/9) : (3.5 x 8/9), (5) 

which lie between the values of Eqs. (|J) and (Q). 



3 Experimental set-up 

The measurements were performed with the ZEUS detector at the HERA ep collider. 
The data, corresponding to an integrated luminosity of 25.0 ± 0.4pb _1 , were collected 
in 1996 and 1997, when HERA operated with a proton-beam energy of 820 GeV and a 
positron-beam energy of 27.5 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [23]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [24], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
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polar- angleQ region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [25] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/^/E for electrons 
and a(E)/E = 0.35/ ' y/E for hadrons (E in GeV). 

A proton- remnant tagger (PRT1) [26] was used to tag events in which the proton dissoci- 
ates. It consists of two layers of scintillation counters at Z = 5.15 m. The pseudorapidity 
range covered by the PRT1 is 4.3 < r\ < 5.8. 

The photoproduction tagger (PT) [10] is a small electromagnetic calorimeter located close 
to the beam-pipe at Z = —44 m. It detects positrons of energy between 21 and 26 GeV 
scattered under very small angles (less than a few mrads). The positron measured in the 
PT is used to tag photoproduction events with a photon-beam energy in the 1.5-6.5 GeV 
range. 

The luminosity is determined from the rate of the Bethe-Heitler bremsstrahlung process 
e + p — > e + ^fp, where the high-energy photon is detected in a lead-scintillator calorimeter 
(LUMI) located at Z = -107 m in the HERA tunnel [27]. 

4 Kinematics and decay distributions 

Figure |l](d) shows a diagram of the reaction 

e(k)p(P) -> e(k')V(v)Y(P'), 

where V is a p°,<f), or J/ip meson and k, k', P, P' and v are the four-momenta of the 
incident positron, scattered positron, incident proton, diffractively produced state Y and 
vector meson, respectively. The kinematic variables used to describe proton-dissociative 
vector-meson production are: 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 
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• Q 2 = —q 2 = —(k — k') 2 , the negative of the squared four-momentum of the exchanged 
photon; 

• y = (P ■ q)/(P ■ k), the fraction of the positron energy transferred to the photon in the 
rest frame of the initial-state proton; 

• W 2 = (q + P) 2 , the squared centre-of-mass energy of the photon-proton system; 

• t = (P — P') 2 = (v — q) 2 , the squared four-momentum exchanged at the hadronic 

vertex; 

• x, the fraction of the proton's momentum carried by the scattered parton. Assuming 
that the exchanged object in the photon-proton scattering couples to a single massless 
parton in the proton 

X = 2P(P~-P) = - tim ~ M * ~ ^ (6) 

where M p is the mass of the proton and M 2 = {P') 2 is the squared mass of the 
diffractively produced state, Y. 

The angles used are shown in Fig. |2| and are: 

• <£>, the angle, in the photon-proton centre-of-mass frame, between the positron-scattering 
plane and the ^-production plane. The latter contains the momentum vectors of the 
virtual photon and the vector meson; 

• 9h and <fh, the polar and azimuthal angles of the positively charged decay particle in 
the V helicity frame. The helicity frame is the rest frame of the meson with the meson 
direction in the photon-proton centre-of-mass frame taken as the quantisation axis. 
The polar angle, 8h, is defined as the angle between the direction of the positively 
charged decay particle and the quantisation axis. The azimuthal angle, (fih, is the 
angle between the decay plane and the ^-production plane. 

Only the momentum vectors of the meson decay products were measured. The momen- 
tum of the scattered positron (and hence Q 2 ), the angle $, and the mass My were not 
measured. 



In the tagged photoproduction events, Q 2 ranges from Qmin ~ 10~ 9 GeV 2 up to Q^ax ~ 
0.02 GeV 2 , with a median Q 2 of approximately 7 x 10~ 6 GeV 2 [10]. The value of Q 2 was 
neglected in the reconstruction of the other kinematic variables. 

The variables W and t can be expressed as [10] 

W 2 w 2E P (E - p z ) and t w -jft,, 
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where E p is the incoming proton energy and E is the energy, pz is the longitudinal 
momentum and pr is the transverse momentum of the produced vector meson in the 
laboratory frame. 

Since $ is not measured, the three-dimensional angular distribution was averaged over 
this angle, so that the ±1 photon-helicity states are not distinguished. The normalised 
two-dimensional angular distribution can then be written in terms of spin density matrix 
elements as 

1 d 2 <r 



=F 



V2Re[r°g] sin 29 h cos cp h T rf_ x sin 2 6 h cos 2<p h } , (7) 



where the upper (lower) signs are for the V decay to spin-0 (spin-1/2) particles. Integra- 
tion over Oh or (p^ produces the one- dimensional distributions 

ocl T r 04 ±(3r ° 4 -l)cos 2 ^ (8) 



dcos^ 
and 

da 



oc 1 =F r^_i cos2tp h . (9) 



dip h 



Since Q 2 is very small for this study, the contribution from photons in the helicity-0 state 
is negligible. In this case, the matrix elements are related to the helicity amplitudes, 
T\ v \ 1 , as follows [28]: 



„04 



3n r p 4 ^ 1 (T n T *) + (T„iiT *) M 1 gugii) + (T-nT*) 



00 ati2 i T^2 i T^2 ' 10 9 rp2 , _i_ T^2 ' 1-1 ~ o T^2 _i_ T^2 _i_ 

01 11 — 11 z 01 11 — 11 z 01 11 — 11 



(10) 



where Tn is the helicity non-flip amplitude, T i is the helicity single-flip amplitude and 
T_ii is the helicity double-flip amplitude. If SCHC holds, only the Tn amplitude is 
non-zero and therefore r[] 4 , and should all be zero. 



5 Event selection 

Events were selected online with a three-level trigger system [23]. The selection of the 
reaction ep — > eVY at large values of —t required: 

• at the first trigger level: 

— a minimum energy deposit of 1 GeV in the PT; 
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— at least one track candidate in the CTD. 

In addition, the following cuts were applied to further reduce the backgrounds from 
random coincidences between some activity in the CTD and a detected positron in 
the PT: 

— an upper limit of 1 GeV on the energy deposited in the LUMI photon detector. This 
requirement suppresses bremsstrahlung events and rejects hard QED radiation; 

— at least 0.5 GeV of energy deposit in the EMC section of the RCAL. This cut 
rejects proton-gas events. However, it also reduces the acceptance for ep — > eVY 
events and therefore was used only when necessary; 

• at the second trigger level: 

— CAL timing consistent with an ep collision; 

— not more than three tracks in the CTD; 

• at the third trigger level: 

— exactly two tracks in the CTD pointing to a common vertex within —60 < Z < 
60 cm; 

— the transverse momentum of at least one track candidate greater than 0.8 GeV. 
This cut efficiently selects events with large values of —t. 

After the event reconstruction, the following requirements were imposed: 

• an interaction vertex within —40 < Z < 40 cm and a transverse distance within 0.7 cm 
of the nominal interaction point (IP); 

• exactly two tracks with opposite charges, each associated with the reconstructed event 
vertex, and each with pseudorapidity \rj\ < 2.1 and transverse momentum greater than 
150 MeV for the p°, 400 MeV for the 0, and 500 MeV for the J/ip candidates; 

• energy deposits in the CAL (excluding FCAL cells with t] > 2.1), not associated with 
tracks, smaller than 250 MeV. This cut rejects events with an additional particle not 
associated with either the vector meson or the proton-dissociative system. 

In addition, -t > 1.1 GeV 2 and 80 < W < 120 GeV were required to select a kinematic 
region of high acceptance. The final number of events, after all selection cuts, was 21414 
for 0.55 < M wn < 1.3 GeV, 2407 for 0.99 < M KK < 1.06 GeV and 214 for 2.7 < M u < 
3.2 GeV. 

6 Monte Carlo simulation and acceptance calculation 
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6.1 Monte Carlo simulation 

The EPSOFT [29] Monte Carlo (MC) generator was used for the simulation of the reaction 
7P — > VY. The particular version [30] used here had an improved simulation of the final- 
state particle multiplicity. The DIFFVM [31] MC generator was used for systematic 
checks. In these generators, 7p interactions are simulated assuming the exchange of 
a colourless object. The cross-section d?a lp ^y Y / (dtdM y ) was parameterised using an 
exponential t distribution and a 1/M-p- distribution; SCHC was also assumed. Because 
the MC samples thus generated did not give a sufficiently good description of the data, 
they were reweighted as described in the next section. The only difference between the 
EPSOFT and DIFFVM generators that is relevant for this analysis is the fragmentation 
scheme of system Y: 

• in EPSOFT, the exchanged object is assumed to couple to the whole proton (Fig. |l|(a)), 
which subsequently fragments into a state Y . The particle multiplicities and the 
transverse momenta of the hadrons in Y are simulated using parameterisations of 
pp data [32], while the longitudinal momenta are generated with a uniform rapidity 
distribution; 

• in DIFFVM, the exchanged object couples to a single quark (Fig. |l](b)) in the proton. 
The fraction of the proton's momentum carried by the struck quark is given by Eq. @. 
The struck quark and the diquark remnant are colour-connected and are fragmented 
into the final-state Y by the JETSET7.3 [33] program. 

As pointed out in Section |2.2| , pQCD models predict that the process 7p — > VY is 
dominated by the reaction involving a gluon (Fig. |](c)) from the proton. This process 
was implemented as an option in DIFFVM. In the following sections, the DIFFVM sample 
where this option was selected is called DIFFVMg, while the DIFFVM sample in which 
the proton emits a quark is called DIFFVMq. 

The response of the ZEUS detector to generated particles was simulated using a pro- 
gram based on GEANT3.13 [34]. The generated events were processed through the same 
selection and reconstruction procedures as the data. 

6.2 Modelling of the dissociative-system Y 

The FCAL energy distributions are sensitive to the fragmentation schemes of the system 
Y as well as to the shape of the generated My distribution. The latter was reweighted 
with a function (1/M y )^*) -1 to have the effective dependence 




(11) 
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The exponent j3(t) was chosen so as to reproduce the measured energy distributions in 
the FCAL for the p° data sample: (3{t) = 0.93e an< for EPSOFT and (3{t) = 0.67e" ait 
for DIFFVMg. These parameterisations are also valid for the and J/ip mesons. The 
comparisons between the data and the MC simulations are displayed in Fig. |3| for several 
t ranges and for the three mesons. The agreement is satisfactory for DIFFVMg and 
EPSOFT but not for DIFFVMq (shown only for the p° in Fig. ||]). This observation is 
consistent with the pQCD expectation that the photoproduction of vector mesons with 
proton dissociation at low x is dominated by the reaction in which a gluon from the 
proton couples to the vector meson. The DIFFVMg samples were used for systematic 
checks since the results do not change significantly when the mixture of DIFFVMq and 
DIFFVMg expected in Eq. ([|) was used. 



6.3 Acceptance corrections 

The overall acceptance is the product of the PT acceptance and that of the main detec- 
tor. The geometric acceptance of the PT was evaluated with a program that uses the 
HERA beam-transport matrices to track the positron through the HERA beamline. The 
simulation was tuned so as to reproduce the tagging efficiency of the PT for Bethe-Heitler 
events, ep — > e^p, triggered by a photon measured in the LUMI photon detector. This 
procedure was described previously [10]. The variations of the position and tilt of the 
positron beam at the IP observed during data taking [35] were also simulated. 

The efficiency for photoproduction events was determined using the geometric acceptance 
of the PT and events were generated according to the equivalent-photon approxima- 
tion [36]. The photoproduction tagging efficiency, calculated as a function of the positron 
energy, E e i, is shown in Fig. £|. The range 80 < W < 120 GeV was chosen so that the 
systematic uncertainty of the tagging efficiency, shown by the shaded area, does not dom- 
inate the total systematic uncertainty. For the kinematic range used in this analysis, the 
cross-section-weighted PT acceptance averages 80 ± 6% Q. 

The decay-angle distributions at the generator level were iteratively reweighted according 
to Eq. (0) using the spin density matrix elements found in this study (see Section |9.6| ) 



The spin density matrix elements for the p° and <fi mesons were assumed to be linear 
functions of t; SCHC was assumed for the J/ip meson. The generated t distributions were 
reweighted with a function e p ^\ where P(t) is a polynomial, to reproduce those observed 
in the data. The resulting simulation agrees reasonably well with the data, as shown 
in Fig. |5|. 



2 For these data, the PT was shifted closer to the positron beam, yielding a higher tagging efficiency 
than in the previous ZEUS measurement [10]. 
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Acceptance corrections evaluated on the basis of the reweighted EPSOFT sample, which 
gives the best description of the data, were calculated separately for each bin of any given 
variable. 

On average, 6% of ep — > eVY events were rejected by the trigger due to accidental coin- 
cidence with bremsstrahlung photons in the LUMI calorimeter. This effect was corrected 
for in the analysis. 

7 Backgrounds 

The dominant background sources, common to all channels, came from the double- 
dissociative, 7p — > XY, and non-diffractive, 7p — > anything, processes. Additional back- 
grounds contributing to individual channels are non-resonant 7r + 7r~ production for the p° 
sample, p° production for the <f> sample, and inelastic Bethe-Heitler production of e + e~ 
and pairs for the J /if) sample. These backgrounds were statistically subtracted 

using the fits to the invariant-mass distributions described in Section \).2\ 

The contribution of elastic vector-meson production (jp — > Vp) to the data was estimated 
from the fraction of events, F, with a signal in one of the PRT1 counters above a threshold 
corresponding to the signal of a minimum-ionising particle. For elastic vector-meson pro- 
duction, such energy deposits are absent. A comparison of the EPSOFT prediction with 
the measured values of F shows that the contribution from the elastic process decreases 
from about 0.2 at -t = 1.2 GeV 2 to about 0.05 at -t = 3.0 GeV 2 . The precision with 
which this background is known is 0.01 for the p°, 0.03 for the and 0.1 for the J /if) 
meson. This background was statistically subtracted bin by bin. For — t > 3 GeV 2 , the 
contamination from the elastic reaction was consistent with zero and was neglected. 

The contributions from the u) and mesons (decaying to 7r + 7r~7r°) in the p°-candidate 
sample were estimated using a MC simulation and the cross-section ratio given by the 
SU(4) prediction of 1/9 for u)/p° and the results of this analysis for 0/p° (see Section [9.4J) . 
This background was less than 3%. The contribution from the if)' meson (decaying to 
J/if>+ neutrals) in the J /if) candidate sample was less than 3%. These backgrounds 
contribute to the overall normalisation uncertainty. 

8 Systematic uncertainties 

The overall systematic uncertainties in the cross sections and density matrix elements in 
each bin were obtained by summing in quadrature the uncertainties listed below: 

• event selection: 
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— the selection cuts, described in Section |], were varied within the resolutions of the 
cut variables; 

— a different procedure to associate CTD tracks and CAL objects was used. 

The effect on the cross sections was less than ±10%, and less than ±0.03 in the spin 
density matrix elements; 

• signal extraction: 

— the parameterisations of the signal shape were changed from the functional forms 
to the MC expectations; 

— the ranges of the M niX) M KK and Mu mass distributions used in the fit procedure 
described in Section fj.2| were varied by their resolutions. 

This resulted in changes of less than ±10% in the cross sections and less than ±0.01 
in the spin density matrix elements; 

• modelling of the dissociative-system Y: 

— P(t) in Eq. (0) was changed to /?(*) = (0.93 ± 0.07)e (ail±a06) * for the p°, (3{t) = 
(0.93±0.17)e( ai9±a17 )* for the and (3(t) = 0.73 ±0.37 for the J/ip meson. These 
ranges were chosen to maintain a satisfactory agreement between the measured 
FCAL energy distributions and the EPSOFT simulation; 

— the fragmentation scheme of the system Y was changed by using DIFFVMg instead 
of EPSOFT for the acceptance corrections; 

— M Y -M%-t was used instead of M Y [37] in Eq. QTI]); 

— the 1/My dependence was reweighted to agree with the pQCD prediction of 
Eq. © Q. 

The resulting uncertainty in the cross sections, coming mainly from the (3(t) variation, 
was less than ±20%. All checks in this class change the generated My distribution. 
Therefore, to avoid double counting, the maximum deviation from the nominal value 
was taken as the total systematic uncertainty of this class of systematic check. Since 
the MC simulation assumes factorisation of the proton and meson vertices, these 
systematic checks result in no changes to the spin density matrix elements; 

• reweighting of the other MC distributions in EPSOFT: 

— the 9h and <fh distributions were reweighted according to the SCHC predictions 
for the p° and mesons. For the J/ip, the values = 0.13 or r^_ 1 = ±0.2 were 
used, corresponding to one standard deviation from the SCHC values, as discussed 
in Section |9~6l ; 

3 The GRV98 [19] parton densities were used, but the results do not change when other parameterisations 
are selected. 
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— the t and W distributions were reweighted by a power-law function of t and W, 
respectively, whilst maintaining satisfactory agreement between the data and the 
MC distributions. 

These modifications to the analysis procedure lead to changes in the cross section of 
less than 7% and less than ±0.01 in the spin density matrix elements; 

• calorimeter trigger efficiency: 

— the EMC RCAL trigger efficiency was varied within its uncertainty, as determined 
from subsamples of events triggered without the requirement of an energy deposit 
in the EMC section of the RCAL. 

This resulted in changes to the cross sections of less than ±5% and variations in the 
spin density matrix elements of up to ±0.01; 

• normalisation uncertainty: 

— PT acceptance; the systematic uncertainty related to the PT acceptance was eval- 
uated by varying the inputs to the MC simulation. The energy scale for photons 
measured in the LUMI detector was changed by ±50MeV [27], the horizontal 
position of the PT by ±0.3mm, and the average position of the IP by ±0.1mm. 
These variations correspond to the systematic uncertainties, as calculated from 
the Bethe-Heitler event sample, discussed in Section |6.3| . The effect on the cross 
sections was ±6%; 

— tracking-trigger efficiency; the uncertainty on the CTD trigger efficiency was esti- 
mated from an independently triggered data sample and corresponds to an uncer- 
tainty on the cross section of ±4%; 

— PRT1 acceptance; the systematic uncertainty related to the PRT1 acceptance was 
estimated from the difference between the fraction of events with a signal in the 
PRT1 in the data and in the EPSOFT simulation in the region of large —t, where 
the elastic contribution is expected to be negligible. This resulted in a systematic 
uncertainty of ±3% in the cross sections. This value was also used in the low —t 
region; 

— other uncertainties; backgrounds from the uj and </> mesons in the p° sample and 
from the ip' meson in the J/ip sample (—3%), as discussed in Section |7] ; QED 
radiation (±2%) [10]; measurement of the integrated luminosity (±1.8%) [27]; 

The overall normalisation uncertainty of ±10% was the sum in quadrature of the above 
uncertainties. This uncertainty does not affect the spin density matrix elements. 
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9 Results 



Although My (and hence x) was not measured, the trigger and offline selection restricted 
M Y to be < 25 GeV, which corresponds to 0.005 < x < 1 at -t = 2 GeV 2 and 0.03 < x < 1 
at —t = 10 GeV 2 . To facilitate the comparison of the data with the pQCD predictions, 
the cross sections presented in the following sections were evaluated in the fixed x range 
0.01 < x < 1, which approximately corresponds to the x range covered by the data at the 
average t of this analysis. 



9.1 Cross-section determination 

The differential cross-section da 7P ^yy/d£ for the photoproduction process 7p — > VY was 
obtained from the cross section measured for the reaction ep —>■ eVY in each bin of the 
variable £(= t, ip, cos Oh) as 

d(X 7 p^yy _ da ep ^ e VY 1 _ Ny 

d£ ~ d£ % ~ zvF^V^V 

where $ 7 is the effective photon flux [36], Ny is the background-corrected and acceptance- 
corrected number of vector mesons in the range A£, extracted using the fits described 
in Section |9.2| , C is the integrated luminosity and By is the branching ratio of the vector- 
meson decay channel. The cross sections for p° production were integrated over the range 
2M V < M n7T <M P + 5T P . 



9.2 Mass distributions 

The 7T7T, KK, and 11 mass spectra for two representative t ranges, together with the results 
of the fits discussed below, are shown in Fig. || 

The acceptance-corrected mr mass spectra for the p°-candidate sample were fitted in the 
range 0.55 < M n7T < 1.3 GeV using the Soding parameterisation [38]: 



da 



A 2 



y/M^MpoTpO 



M^-M 2 p0 + iM p0 T p o 



+ B/Afp 



+ fps(M v 



(12) 



where B is the non-resonant amplitude (taken to be constant and real), A p o is the normal- 
isation of the resonant amplitude, M p o is the p° mass and r p o is the momentum-dependent 
p° width [39]. The additional term f PS = A PS (1 + 1.5M n7T ) [10] was used to account for 
residual background from the non-diffractive ^p — > anything and the diffractive ^p — > XY 
processes. The ^/ndf for all the fits is satisfactory. The ratio \B/A p o\ decreases from 
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0.13 ± 0.02GeV~ 1/2 at -t = 1.3 GeV 2 to 0.05 ± 0.03GeV" 1/2 at -t = 2.3 GeV 2 and is 
consistent with zero above —t = 2.5 GeV 2 . The background contribution in the range 
0.6 < M n7T < 1.3 GeV, given by the integral of the function fp$, increases from about 
8 ± 2% at -t = 1.3 GeV 2 to 17 ± 4% at -t = 2.7 GeV 2 and is consistent with 15% above 
-t = 3.2 GeV 2 . 

The detector efficiency does not vary over the region of the <fi and J / ip peaks. The <ft and 
J/tp signals were, therefore, extracted from the uncorrected mass spectra. The KK mass 
spectra for the sample of 4> candidates were fitted in the range 0.99 < Mkk < 1-15 GeV 
with the sum of a Breit-Wigner function convoluted with a Gaussian resolution function 
for the signal, and the function (Mkk—%Mk) s for the background. The x 2 / n df for all the 
fits is satisfactory. The background contribution in the range 0.99 < M KK < 1.06 GeV 
decreases from 25 ± 2% at -t = 1.3 GeV 2 to 3 ± 2% above -t = 4 GeV 2 . 

The 11 mass spectra for the J/ip candidate sample in the dilepton- invariant-mass range 
2.6 < Mu < 3.4 GeV were fitted with the sum of a Gaussian (for the muon decay channel), 
a bremsstrahlung function convoluted with a Gaussian (for the electron decay channel), 
and a background term linear in Mu. No muon or electron identification was performed; 
the relative contributions of muon and electron pairs were taken from the MC simulation. 
The electron mass0 was used in calculating Mu. The x 2 /ndf for all the fits is satisfactory. 
The background contribution in the range 2.7 < Mu < 3.2 GeV is 25 ±2%, approximately 
independent of the t range under study. 

Since the backgrounds depend on t, W, cos 9h and <ph, each of the results presented in the 
following sections was obtained by repeating the fits to the invariant mass distributions 
in each kinematic bin. 

9.3 Measurement of da 7P -+vY / 'dt 

Tables [I] - [3] and Fig. [7] show the differential cross sections for the p°, and J/ip mesons 
obtained from the present data together with results obtained in the ZEUS 1995 [10] 
study of the lower —t region. In that analysis, the cross section was determined for 
M^/W 2 < 0.01. In the kinematic region -t « 1 GeV 2 and W ps 100 GeV, this limit 
is numerically equivalent to that used in the present study, so that the results may be 
directly compared. The two measurements are consistent in the region of overlap. 

The cross sections are well described by a power-law dependence, as expected for a hard 
production mechanism [9]. A fit to the present data with the function A ■ (— t)~ n gives 
n = 3.21 ± 0.04(stat.)±0.15(syst.) for the p°, n = 2.7 ± 0.1(stat.)±0.2(syst.) for the 0, 
and n = 1.7 ± 0.2(stat.)±0.3(syst.) for the J/ip meson. 

4 The results do not change when the muon mass is used. 
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The measured values of n for the p° and (J) mesons are significantly smaller than the lower 
limit n ~ 3.8 predicted by Ivanov et al. [17]. The value of n for the J/ip meson agrees 
with the Bartels et al. [16] BFKL prediction of 1.8. The BFKL prediction of Forshaw 
and Poludniowski [18], for the most natural choice of a s = 0.2 and Wq = M v , gives 
X 2 /ndf = 0.5 for all three mesons and is shown in Fig. || These measurements establish 
that perturbative QCD calculations are able to describe the large —t regime. The dip at 
—t ~ M v predicted by the two-gluon exchange model of Forshaw and Ryskin [13] is not 
seen [18] in the data for light mesons, indicating that this model is not applicable in the 
— t ~ My region. 

Note that n depends slightly on the x range used in the cross-section definition both in 
the experiment and in the theory. Hence the values quoted in this paper apply only to 
the cross section integrated over x > 0.01. 



9.4 Ratios of cross sections 

The cross-section ratios 

r * = dt 1 dt 

and 

_ da^p^j/^y /d<J 7 p^pOy 
rj/ + = dt ' dt 

are shown in Tables |] and |5] and in Fig. |9|, together with the values given in Eqs. ([3]), (|J) 
and (||). A clear increase of and rj/^ with increasing —t is observed. 

The dependence of and rj/^ on t suggests that the production mechanism may become 
flavour independent at larger —t values than covered in this analysis. In the framework 
of pQCD, this observation supports the hypothesis that the quark mass appearing in the 
'-/V coupling should be the constituent quark mass, since, in this case, the asymptotic 
value of should be reached at —t ^> M^ , M|. If the quark mass is interpreted as the 
current quark mass, is expected to approach its asymptotic value at much smaller —t 
values. 

The cross-section ratios given in Eq. expected for the production of mesons in the 
helicity-0 state by transverse photons at large t [15,17], are the same as those expected 
from pQCD for the elastic production of mesons in the helicity-0 state by longitudinal 
photons at large Q 2 [8,40]. The cross-section ratios for production of mesons in the 
helicity-±l states are expected to be larger than those for production of mesons in the 
helicity-0 state. Figure |1(] compares the t and Q 2 dependences of the ratios for the 
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proton-dissociative and for the elastic process, respectively; lies systematically above 
cr-y'p^tjyp/cr^p^pOp [41,42] but is compatible within the experimental uncertainties, while 
rj/ip for —t > 3 GeV 2 is about a factor of five larger than a^p^j/^p/ a^*p^ p o p for Q 2 > 
3 GeV 2 as measured by ZEUS [43] and HI [44, 45] f\ The helicity analysis, described in 
Section |9.6| , shows that the production of mesons in the helicity-±l states dominates in 
the large —t region. 



9.5 The slope of the Pomeron trajectory 

The cross section given by Eq. ([I]), integrated over My, can be expressed as 

where F(t) is a function of t only. Assuming a linear form for the Pomeron trajectory 
implies 

da Y VY = F (t)W Aa ^ +Wt - A = F(t)G(W)W 4a '\ (13) 

(-16 

where G(W) is a function of W only. 

The absolute W dependence, and hence a(t), cannot be measured using the present data 
because of the large anti-correlated systematic uncertainties of the PT acceptance as a 
function of W. However, these uncertainties do not affect the measurement of the relative 
changes of the W dependence as a function of t, which are sensitive to a'. The slope of the 
Pomeron trajectory for the p° meson was obtained by fitting the form given by Eq. ( |l3|) to 
the measured cross sections using three W times six t bins. A two-dimensional fit gave a' = 
-0.02 ± 0.05(stat.)+J;Jg(syst.) GeV" 2 with x 2 /ndf = 4.6/9. A similar fit for the meson, 
performed in three W times five t bins, gave a' = -0.06 ± 0.12(stat.)+S;Sg(syst.) GeV~ 2 
with x 2 / n df = 3.4/7. There were insufficient data to perform a similar analysis on the 
J/ip sample. 

The consistency of the a' values with zero implies that the W dependence does not vary 
with t, and vice versa. This conclusion is manifest when looking at the W dependence of 
the ratio of differential cross sections 

dcr{W) 

dt l*=*o 



shown in Figs. |IT] and |T2|. Values of — 1 = 1.3 GeV were chosen for the p° and —t^ 
1.5 GeV 2 for the 6 case. 



The J/i(j data from HI have been scaled from W — 90 GeV to W = 75 GeV and the p° data calculated 
at the same Q 2 values using a fit of the form (Q 2 + M 2 a )~ n [46] 
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Figure 13 compares the values of a' measured in this analysis with those obtained at 
lower —t. For both p° and <p mesons, the values of a' for —t > 1.3 GeV 2 are smaller 
than the value a' = 0.25 GeV -2 characteristic of soft hadronic processes at —t < 0.5 GeV 2 
and also below the values a' = 0.125 ± 0.038 GeV" 2 , a' = 0.158 ± 0.028 GeV" 2 and 
a' = 0.113 ± 0.018 GeV" 2 obtained for -t < 1.5 GeV 2 , from fits [10, 11] to the elastic 
photoproduction of p°, <p an d J ftp mesons, respectively. This observation indicates that 
a' decreases with increasing —t. The small values of a' obtained in this analysis are 
consistent with the flatness of the Pomeron trajectory observed by UA8 [47] at large —t 
and with the expectation of the anomalous Regge exchange model by Kochelev et al. [12]. 
They are also consistent with the expectation of a' = characteristic of the two-gluon- 
exchange models [13,15,17,18] and the expectation of a' < 0.1 GeV" 2 of the BFKL model 
of Bartels et al. [16]. 



9.6 Decay angular distributions 

The r°Q and r^z 1 spin density matrix elements for the p°, <p, and J/ip were obtained by 
fitting Eqs. (§) and (|) to the background-subtracted and acceptance-corrected data in 
several t ranges. The t ranges and the results of the fits are given in Tables |6|, |^ and 
H As an example, Fig. [14| shows the cos 9h and tfh distributions with the results of the 
fits in the highest —t bin used in the helicity analysis for each of the three mesons. The 
flat (fih distribution shown by the dashed line disagrees (x 2 /ndf = 42/7) with the p° 
data, indicating a violation of SCHC. The parameter Refr^] for the p° and </> mesons was 
determined using a two-dimensional fit of Eq. (^) in three and two t intervals, respectively. 
The results of the fits are summarised in Tables |6] and 0. The above results were also 
confirmed by a decay-angle analysis [48] performed in the transversity frame 0. There 
were insufficient data to perform a two-dimensional fit for the J/ip meson. 

The fitted values of Tqq, Refr^], and for the p° meson are displayed in Fig. |l^. The 
results are consistent with the ZEUS 1995 [10] measurements in the overlap region. The 
small values of indicate that the probability to produce p° mesons in the helicity-0 



state from a photon in the helicity-±l states is 4 ± 5% at —t = 3.35 GeV and 9 ± 
at —t = 5.67 GeV 2 . However, the non-zero value of Re[r°Q] indicates a helicity single-flip 
contribution at the level of a few percent. The finite negative values of show clear 
evidence for a helicity-double-flip contribution. 

The measurements for the meson, shown in Fig. [16|, display similar features to those 
for the p° meson, although with a smaller statistical significance. The values of r°Q and 



6 The spin-quantisation axis in the transversity frame is chosen along the normal to the meson production 
frame. 
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rjij for the J/ip meson, shown in Fig. 17, are consistent with the SCHC hypothesis. 

Although the BFKL calculation of Forshaw and Poludniowski [18] successfully describes 
the cross sections for all three mesons, it is not in accord with the observation of SCHC 
violation. This calculation forces SCHC by using a delta-function for the meson wave- 
function. 

If a more appropriate meson wave-function for the light mesons is used, the helicity of the 
final-state meson will depend on the modelling of the photon fluctuation into a qq pair. 
Within perturbation theory, a photon can only split into a qq pair in a chiral-even spin 
configuration. If, in addition, the quark appearing in the jV coupling is interpreted as a 
current quark, as assumed by Ginzburg and Ivanov [15], real photons will produce light 
mesons only in the helicity-0 state because of the chiral nature of the perturbative coupling 
in the massless limit. Not only does this model not agree with the cross section [10], 
but it is also ruled out by the helicity data. Ivanov et al. [17] have proposed that a 
large contribution of light mesons produced in the helicity-±l states could arise from 
qq fluctuations in a chiral-odd spin configuration. This model fails to describe the cross 
section and does not quantitatively describe the helicity structure. 



10 Summary 

Diffractive photoproduction of p°, <p and J/ip mesons with proton dissociation has been 
measured in the ZEUS detector at HERA using data at W ~ 100 GeV for —t values up to 
12GeV 2 . Expressing the differential cross section as a power-law, da^^yy/dt oc (— t)~ n , 
yields a value for n that decreases with increasing mass of the meson. The resulting expo- 
nents are: 3.21±0.04(stat.)±0.15(syst.) for p° production; 2.7 ± 0.1(stat.)±0.2(syst.) for 
(f> production; and 1.7 ± 0.2(stat.)±0.3(syst.) for J/ip production. The BFKL calculation 
of Bartels et al. [16] agrees with the J/ip data. An extension of this model by Forshaw 
and Poludniowski [18] is able to describe the data for all three mesons. 

The cross-section ratios for the <fi and J/ip with respect to the p as a function of t in- 
crease with —t. In the context of pQCD, this observation suggests that the production 
mechanism may become flavour independent at larger —t values than those covered in 
this analysis. 

The decay-angle analysis for the p° and <fi mesons indicates that the SCHC hypothesis 
does not hold, since both single and double helicity-flip contributions are observed. No 
available pQCD calculation is able to describe this result quantitatively. 

The effective slope of the Pomeron trajectory, a', was determined by studying the relative 
changes of the W dependence of the cross-section as a function of t for the p° and <fi 
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mesons. The slopes are a' = —0.02 ± 0.05(stat.)^Q og (syst.) GeV~ 2 for the p° meson and 
a' = — 0.06±0.12(stat.)^QQg(syst.) GeV~ 2 for the <t> meson. These values are lower than the 
value characteristic of low —t hadronic scattering (a' = 0.25 GeV -2 ) and lower than those 
obtained in elastic vector-meson photoproduction at HERA at lower —t [10,11]. They 
are in agreement with pQCD expectations [16, 17] and are consistent with the flatness of 
the Pomeron trajectory at large —t observed by UA8 [47]. These results imply that large 
values of —t can provide a suitable hard scale for perturbative QCD calculations. 
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ZEUS 1996-97 1P -> p°Y 


—t range 
( GeV 2 ) 


(-*) 

\ / 

( GeV 2 ) 


da/d\t\ 
(rib/ GeV 2 ) 


1.20-1.50 


1.33 


163 ± 6 + \t + \ 2 

— Li) — 1Z 


1.50-2.00 


1.71 


80.1 ± 2.5 ±|:| tl° 5 


2.00-2.50 


2.22 


36.9 ± 1.4 tli t° 2 i 


2.50-3.00 


2.72 


17.5 ± 0.9 t\i t?; 7 


3.00-3.50 


3.23 


10.7 ± 0.7 t° 8 7 t°{ 3 2 


3.50-4.00 


3.73 


6.98 ± 0.59 ±£g 


4.00-5.00 


4.43 


3.39 ± 0.27 j£™ + -°o% 


5.00-6.85 


5.75 


141 i 15 +0 ' 35 + - 20 

l.ftl n= u.io _ 049 _ 025 


6.85-10.0 


8.04 


40 + 07 +0 ' 17 + - 08 

U.^U 3= U.U/ _ 022 _ .13 



Table 1: T/ie differential cross section, dcr 7P ^ p oy/d|t| ; /on > 0.01 and 80 < W < 
120 (7eU. The first uncertainty is the statistical, the second the systematic and the 
last is due to the modelling of the proton-dissociation process. The normalisation 
uncertainty of 10% is not included. 



ZEUS 1996-97 ~fp -> 0Y 


—t range 
( GeV 2 ) 


(-*) 
( GeV 2 ) 


d*/d|t|) 
(nb/ GeV 2 ) 


1.20-1.55 


1.36 


19.2 ± 1.1 tli j£? 


1.55-1.90 


1.71 


12.4 ± 0.7 t\i t°oi 


1.90-2.30 


2.08 


7 40 i n +0.68 +0.30 
' 31 u '°° -0.72 -0.41 


2.30-2.70 


2.48 


A KQ 4- n AC) +0-49 +0.17 
^.OO ^ u -^ u -0.43 -0.08 


2.70-3.20 


2.93 


2 ^ + 26 +0 - 30 +0 - 15 

Z.OO HZ U.ZU _Q 2 g _o.08 


3.20-3.80 


3.47 


1.55 ± 0.20 ±ss ±S:S 


3.80-4.80 


4.24 


o.89 ± 0.18 ±8:S 


4.80-6.50 


5.53 


0.36 ± 0.08 US 



Table 2: TTie differential cross section, do" 7P ^y/d|t|, /on > 0.01 and 80 < W < 
120 GeU. TTie first uncertainty is the statistical, the second the systematic and the 
last is due to the modelling of the proton-dissociation process. The normalisation 
uncertainty of 10% is not included. 
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ZEUS 1996-97 7 p -> J/^Y 


—t range 


(-*> 


dcr/d|t| 


( GeV 2 ) 


( GeV 2 ) 


(nb/ GeV 2 ) 


1.20-1.50 


1.34 


16.0 ± 2.4 tH 


1.50-2.00 


1.74 


9.5 ± 1.5 t\i t?J 


2.00-2.50 


2.24 


8.3 ± 1.5 t\l ±S5 


2.50-3.25 


2.84 


4 98 ± 9 s ) +0 - 89 +0 - 38 
^.yo in u.jo _ 90 _ 63 


3.25-4.50 


3.81 


2 5 + 06 +0 ' 8 +0 - 4 

Z.O 3= U.U _1.0 -0.4 


4.50-6.50 


5.37 


1 QQ + n +0-80 +0-16 

1.0-3 3= U.^±0 _ Q41 _ 039 



Table 3: TTie differential cross section, da w -,j/^y/d\t\, for x > 0.01 and 
80 < IU < 120 GeU. T/ie first uncertainty is the statistical, the second the sys- 
tematic and the last is due to the modelling of the proton-dissociation process. The 
normalisation uncertainty of 10% is not included. 



ZEUS 1996-97 


—t range 


(-*> 


% (7P-0Y)/^( 7P ^p°y) 


1.20-1.55 


1.35 


0.124 ± 0.008 j™» 


1.55-1.90 


1.71 


0.156 ± 0.011 t° ™ 


1.90-2.30 


2.08 


nivfi + nnK +0.021 +0.011 

U.I/O 3= U.UIO _ .017 -0.010 


2.30-2.70 


2.48 


Oifi7 + ooifi +°- 013 +0- 006 

U.ltW 3= U.U1U _ .013 -0.006 


2.70-3.20 


2.93 


0.185 ± 0.021 ±°Zt Toll 


3.20-3.80 


3.47 


0.189 ± 0.027 t° f 3 °s Toll 


3.80-4.80 


4.23 


223 ± 47 +0 - 057 +0 - 034 

U.ZZO 3= U.U1H -0.042 -0.033 


4.80-6.50 


5.49 


n OOc; 4- nnw +0.055 +0.055 
U.ZZO 3Z U.UOZ _ 0060 -0.052 



Table 4: T/ie raizo 0/ i/ie differential cross sections, da/dt, for <f> to p° proton- 
dissociative photoproduction for x > 0.01 and 80 < W < 120 GeV . The first 
uncertainty is the statistical, the second the systematic and the last is due to the 
modelling of the proton-dissociation process. 
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ZEUS 1996-97 


—t range 


H> 


g( 7 p- J/^Y)/%( 1P ^p°Y) 


1.20-1.50 


1.33 


0.099 ± 0.015 ±™H t° Z 


1.50-2.00 


1.71 


0.119 ± 0.019 $™ ±S 


2.00-2.50 


2.22 


226 + 042 +0 ' 056 +0 - 005 

U.ZZU 3Z U.U^Z _ .037 -0.008 


2.50-3.25 


2.82 


32^ + 63 +0 - 020 +0 - 013 

U.OZO 3= U.UUO _ .059 -0.018 


3.25-4.50 


3.75 


n Ann 4- n 0Q7 +0.078 +0.038 
U.^UJ 3= U.UJI -0.121 -0.046 


4.50-6.50 


5.28 


72 + 2^ +0 - 20 +an 

UJZ 1 U - Z0 -0.22 -0.15 



Table 5: The ratio of the differential cross sections, da/dt, for J/ip to p° proton- 
dissociative photoproduction for x > 0.01 and 80 < W < 120 GeV . The first 
uncertainty is the statistical, the second the systematic and the last is due to the 
modelling of the proton-dissociation process. 



ZEUS 1996-97 7 p -> p°Y 


—t range 
(GeV 2 ) 


H> 

(GeV 2 ) 


r 04 
'00 


Re [r? 4 ] 


r 04 
'1-1 


1.1-1.7 


1.34 


0.022±0.015lg;{£f 


0.054±0.012lHlo 


-0.143±0.014+°;° 2 ? 


1.7-2.9 


2.14 


0.042±0.017lH2i 


0.045±0.013±g;gig 


-0.158±0.017lg;g22 


2.9-4.0 


3.35 


0.037±0.034+°;°| 




-0.129±0.032+°;° 23 


4.0-12 


5.67 


0.090±0.049lg;g^ 




-0.252±0.047lg;g^ 


2.9-12 


4.38 




0.047±0.020lg;g?? 





Table 6: The spin density matrix elements for ^p — > p Y for 80 < W < 120 GeV. 
The first uncertainty is the statistical and the second the systematic. 



ZEUS 1996-97 ~fp -> 0y 


— t range 
(GeV 2 ) 


H> 

(GeV 2 ) 


r 04 
' 00 


Re [r%] 


r 04 
' 1-1 


1.2-1.7 


1.42 


0.080±0.028+°;° 2 ° 


0.053±0.021+°;°^ 


0.008±0.033ig;^ 


1.7-3.0 


2.20 


0.066±0.025inii 




-0.085±0.032l°;° 2 | 


3.0-7.0 


4.03 


-0.020±0.035lg;g?| 




-0.108±0.058ig;g^ 


1.7-7.0 


2.72 




0.057±0.018lg;gJ? 





Table 7: The spin density matrix elements for ^p — » <f>Y for 80 < W < 120 GeV . 
The first uncertainty is the statistical and the second the systematic. 
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ZEUS 1996-97 7 p -> J/^F 


—t range 
(GeV 2 ) 


H> 

(GeV 2 ) 


r 04 
'00 


r 04 
' 1-1 


1.1-1.8 


1.42 


-0.28±0.26lg;jjf 


0.10±0.15l°£ 


1.8-6.5 


3.01 


0.22±0.30lg^? 


-0.11±0.16i°;°4 



Table 8: The spin density matrix elements for ^p — > J/ipY for 80 < W < 
120 GeV . The first uncertainty is the statistical and the second the systematic. 
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p(P) Y(P') ^ 



Figure 1: (a) —> VY scattering in the VDM model. The photon fluctuates 
into a vector meson, which then scatters off the proton and the proton breaks up. 
(b,c) Example of ^p — * VY scattering in pQCD models. The photon fluctuates 
into a qq pair, which then scatters off a single quark (b ) or gluon (c) in the proton 
by the exchange of two gluons or a gluon ladder. The scattered qq pair becomes a 
vector meson. The struck parton and the proton remnant fragment into a system 
of hadrons. (d) Schematic diagram of proton- dissociative vector-meson production 
in ep interactions, ep — > eVY . 



28 



/ / 




Figure 2: Illustration of the angles used to analyse the helicity states of a vector 
meson (for a decay into two particles, V — > q + q~ ). 
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Figure 3: The distributions of the FCAL energy for: (a,b,c) p° ; (d,e,f) <p; (g,h) 
J/ip candidate samples in different —t regions. The symbols represent the data and 
the histograms indicate the results of the simulation using EPSOFT (solid line), 
DIFFVMg (dashed lines) and DIFFVMq (dotted lines shown only for the p°). The 
shaded bands represent the size of the correlated uncertainties due to the modelling 
of the proton dissociation in EPSOFT. 
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Figure 4: The photoproduction tagging acceptance as a function of the energy 
of the scattered positron calculated using simulated events generated according to 
the equivalent-photon approximation in the Q 2 < 0.02 GeV 2 range. The shaded 
band represents the relative systematic uncertainty of the tagging acceptance. The 
dashed lines indicate the kinematic region used in this analysis. 
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Figure 5: Comparison between the data and the EPSOFT MC (histograms) events 
for: (a,b,c) —t; (d,e,f) cos^; (g,h,i) ip^; (j,k,l) W , after all selection cuts (except 
the cut on W in the case of W distributions). The MC is normalised to the data. 
The dashed lines (j,k,l) indicate the cuts on W applied to select the final event 
sample. The systematic uncertainties due to the tagging acceptance are relevant in 
the case of W distributions and are shown as the shaded bands. The three columns 
refer to the p° , <fi, and J/ip samples, as indicated. 
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Figure 6: ('o,^ ?7ie differential cross-section da/dM nn for the p° sample. (c,d) 
The Mkk mass distributions for the <fi sample. (e,f) The Ma mass distributions 
for the J /ip sample. The symbols are the data for representative t ranges and the 
solid curves indicate the result of the fits discussed in the text. The dashed curves 
show the background contributions. 
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Figure 7: TTie differential cross sections der 7P ^yy/d|t| m the range 80 < W < 
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the statistical uncertainty and the outer bars represent the statistical and system- 
atic uncertainties added in quadrature. The shaded bands represent the correlated 
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lated uncertainty of ±10% is not shown. The open symbols correspond to the ZEUS 
1995 results [10]. The lines are the results of the fit to the data with the function 
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Figure 8: The differential cross sections da lp ^vY / d\t\ for: (a) p° , (b) <fi, and 
(c) J/tp photoproduction in the range 80 < W < 120 GeV and x > 0.01. The 
solid curves show a pQCD calculation [18]. The inner bars indicate the statistical 
uncertainty and the outer bars represent the statistical and systematic uncertainties 
added in quadrature. The shaded bands represent the size of the correlated uncer- 
tainties due to the modelling of the hadronic- system Y. An additional correlated 
uncertainty o/±10% is not shown. 
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Figure 9: (a) The ratio of the cross sections da lp ^yy / dt forcf) to p° photoproduc- 
tion. (b) The ratio of the cross sections da lp ^vY / dt for J/ip to p° photoproduction. 
The inset shows this ratio on a linear scale. The inner bars indicate the statistical 
uncertainty, the outer bars represent the statistical and systematic uncertainties 
added in quadrature. The shaded bands represent the size of the correlated un- 
certainties due to the modelling of the dissociative system, Y . The dashed lines 
correspond to the SU(4) predictions, while the dotted and dashed-dotted correspond 
to the pQCD values given by Eqs. (j^) and (^), respectively. The open circles are 
the ZEUS 1995 results [10]. 
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Figure 10: (a) The ratio of the <p to p° cross sections as a function of —t or 
Q 2 . The 4>/p° results as a function of —t for proton-dissociative photoproduction 
from this analysis are shown with solid circles and those from the ZEUS 1995 [10] 
measurement with the solid squares. The shaded bands represent the size of the cor- 
related uncertainties due to the modelling of the dissociative system, Y. Open tri- 
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dashed lines correspond to the SU(4) predictions, while the dotted and dashed-dotted 
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Figure 11: The ratios of p° production cross sections, Tjf/^if |t=t > f or = 1-3 
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Figure 12: The ratios of (p production cross sections / ^\t=t , for —t = 1.5 
GeV 2 , as a function of W in four t intervals. The lines represent the result of the 
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Figure 13: Comparison of results on a' for vector-meson production. The result 
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Figure 14: The normalised background-subtracted and acceptance- corrected cos 8 \ 
and fh distributions for proton-dissociative photoproduction of: (a,b) p°, (c,d) <fi 
and (e,f) J/ip. The symbols represent the data and the solid curves the result of 
the fits with Eqs. (TJ^ and The dashed curves are the SCHC predictions. 
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Figure 15: The fitted values of (a) r^, (b) i?e[r°Q] and (c) r°i x for proton- 
dissociative p° photoproduction as a function of —t. The inner bars indicate the 
statistical uncertainty and the outer bars represent the statistical and systematic 
uncertainties added in quadrature. The open circles correspond to the ZEUS 1995 
results [10]. The SCHC prediction is shown as the dashed line. 
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Figure 16: T7ie /itted values of (a) r^, i?e[r°Q] and (c) r^_ 1 for proton- 
dissociative <p meson photoproduction as a function of —t. The inner bars indicate 
the statistical uncertainty and the outer bars represent the statistical and systematic 
uncertainties added in quadrature. The SCHC prediction is shown as the dashed 
line. 
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Figure 17: The fitted values of (a) Tqq and (b) r°i 1 for proton- dissociative J ftp 
meson photoproduction as a function of —t. The inner bars indicate the statistical 
uncertainty and the outer bars represent the statistical and systematic uncertainties 
added in quadrature. The SCHC prediction is shown as the dashed line. 
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